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The pair localization length L2 of two interacting electrons 
in one-dimensional disordered systems is studied numerically. 
Using two direct approaches, we find L2 oc Lf , where L\ is the 
one-electron localization length and a ~ 1.65. This demon- 
strates the enhancement effect proposed by Shepelyansky, but 
the value of a differs from previous estimates (a = 2) in the 
disorder range considered. We explain this discrepancy using 
a scaling picture recently introduced by Imry and taking into 
account a more accurate distribution than previously assumed 
for the overlap of one-electron wavefunctions. 

PACS numbers: 71.30, 72.15.R 



Very recently, Shepelyansky |l]] considered the prob- 
lem of two interacting electrons in a random potential, 
defined by the Schrodinger equation 



(E - V„ 1<n2 )ip ni! 
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the on-site interaction and Vi is a random potential dis- 
tributed uniformly in the interval [-W, W]. The indices 
n\ and denote the positions of the first and the second 
electron, respectively. Shepelyansky proposed that, as a 
consequence of the interaction U, certain eigenstates ex- 
tend over a range L2 much larger than the one-particle 
localization length L\ oc W~ 2 . A key quantity in the 
derivation of this spectacular result in one dimension 
is the matrix representation IA of the Hubbard inter- 
action in the disorder diagonal basis of localized one- 
electron eigenstates. With i?„^ the amplitude at site n 
of the one-particle eigenstate with energy Ei we have 

Uij,lm — UQijjm with 
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The Qij^im vanish unless all four eigenstates are roughly 
localized within the same box of size L\. Assuming that 
Rn t i oc a n j \fL~\ inside the box, with a n a random num- 
ber of order unity, and neglecting correlations among the 
a n at different sites n, one finds Q w for a typ- 

ical nonvanishing matrix element. In JjJ, this estimate 
was adopted and used to reduce the original problem to 
a certain band matrix model, eventually giving L2 oc L\. 
Later, Imry employed the Thouless scaling block pic- 
ture to reinforce, interpret and generalize this result. The 



key step in this approach involves the pair conductance 
92 = (UQ/S2) 2 , where 82 is the two-particle level spacing 
in a block of size L\ and Q is evaluated between adja- 
cent blocks. Using Q w l/L^ 2 as before, Imry finds that 
g2 ~ 1 on the scale L2 oc L\, in agreement with Shep- 
elyansky. As a second important result both approaches 
predict that the effect does not depend on the sign of U. 

In this letter, we confirm the enhancement effect by 
studying both the original model (0) for finite size sam- 
ples and an infinite "bag model" with medium-range in- 
teraction. However, we find L2 oc with a « 1.65 
instead of a = 2 in both cases. Moreover, the sign of U 
is not entirely irrelevant. We suggest that the small value 
for a is due to a very peculiar distribution of the Qijj m - 
Imry's scaling approach, modified to take into account 
the correct Q-distribution, is shown to be more or less 
consistent with a ~ 1.65 in one dimension. 

Transfer matrix approach for finite size chains. Eq.(^) 
is formally identical to a single-electron Anderson model 
in two dimensions with symmetric (with respect to the 
diagonal n± — 712) on-site disorder. The indices n\ and 
ri2 can be interpreted as cartesian coordinates on a fi- 
nite square lattice with N 2 sites. The iV 2 eigenstates of 
the Hamiltonian are either symmetric or antisymmetric 
with respect to the interchange of m and n.2. The inter- 
action is relevant for the symmetric states only (corre- 
sponding to electrons with opposite spins). For the time 
being we disregard questions of symmetry and proceed in 
analogy to the 2d Anderson model jj[[IJ . Imposing hard- 
wall boundary conditions in the ri2-direction, transfer in 
the ni-direction is described by a product of N matrices 
M p . This amounts to calculating the transport of one 
electron along the Id chain while the other one occupies 
some state (or "channel" ) of the same chain. Note that 
only the energy E of the pair of electrons is fixed. In nu- 
merical computations we are restricted to finite sizes, M p 
being a N x N matrix. This limitation is inessential for 
L2 < N and induces only minor finite-size corrections for 
the transfer matrix eigenvalues otherwise || . We rewrite 
(fll) in the form 
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where u p is a vector with N components (ip P ,i, 
and the transfer matrix M p is defined by: 



M p = 



(E — V p )l — Hi -1 
1 
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Here, Hi is the (transverse) single-particle Hamiltonian 
for the second electron on the disordered chain of TV sites 
with hard-wall boundary conditions, and A p is a diago- 
nal N x N matrix with (A p )a = 6i p . 

Before we numerically evaluate the set of decay lengths 
characterizing this transfer matrix product, we insert a 
few qualitative remarks. With R the orthogonal N x N 
matrix diagonalizing Hi (i.e., Hi = RHdR T , where Hd 
is a diagonal matrix with entries E a ) we transform M p 
according to 



M v = 



' R T 
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In the noninteracting case (U — 0) A4(N) = n^Li-^p 
reduces to N decoupled products of 2 x 2 matrices and 
we deal with a system of N decoupled chains, all de- 
scribed by the same Hamiltonian Hi. On each chain a 
the first electron propagates with energy E — E a while 
the second one occupies an eigenstate of Hi with en- 
ergy E a . The transfer matrix A4 P defines N localization 
lengths £i(E — E a ) related to Hi, eventually becoming 
a continuous density in the thermodynamic limit. In the 
interacting case (U ^ 0) we have to consider in addition 
the N x N matrix UR T A p R for each M p , which gives 
rise to two modifications. First, this matrix couples the 
N chains by a hopping term UR^ p R p p proportional to 
the overlap (at site p) of two eigenstates of Hi with en- 
ergies E a and Ep, respectively. We recover the central 
point of the present problem: the effect of the (local) 
interaction is determined by the spatial overlap of the 
eigenstates of Hi. If these states have random phases, 
the sign of U will be irrelevant in the hopping terms. 
However, the second modification due to U, a shift of V p 
on chain a by UR^ p R p a , will remain sensitive to this 
sign. In short, the interaction couples the Id chains and 
modifies (smoothens or enhances) the fluctuations of V p 
at those sites where the eigenstates of Hi are localized. 

We have numerically studied M(N) — Y\ p= i M p as a 
function of U and W at fixed E. We focus attention 
on the logarithms ji(N) = l/£i(N) of the eigenvalues 
of M T (N)M(N) divided by 2N . One expects this set 
of inverse decay lengths l/£j(iV) to converge to a well- 
defined large- ./V limit We denote by L 2 = £i the 
largest of these lengths to characterize the extension of 
the least localized quantum state of the electron pair. In 
Fig.l we present the five smallest and five largest inverse 
decay lengths 7* (AT) for a given sample of size N = 100 
as a function of U. The curves exhibit large sample- 
dependent fluctuations, but having analyzed many dif- 
ferent samples we can draw the following conclusions: 
(i) the largest lengths are increased by the interaction, 
confirming the effect discovered by Shepelyansky; (ii) as 
expected for local interactions, shorter decay lengths are 
less sensitive to U since the electrons tend to be local- 
ized with small or vanishing overlap; (iii) on average the 



curves are symmetric around U = at the band centre 
E = (due to an exact symmetry of the ensemble of 
Hamiltonians), but this symmetry is broken for E ^ 0, 
and the sign of U appears to be no longer irrelevant; (iv) 
the nearest-neighbour spacing distribution of the l/£j 
shows a crossover from Poisson statistics (for U = 0) to 
Wigner-Dyson statistics at finite U. 

To determine the dependence of Li on W we show in 
Fig. 2 (note the log scale) the ensemble-average L 2 (N) — 
l/<7i (AO) as a function of W for N = 100, E = 0, and 
[7 = 1. One gets L 2 cx W~ 2a with 2a = 3.29, a value 
which differs from the estimate of Shepelyansky [Q and 
Imry 0] (a = 2). Although this qualitatively confirms 
the enhancement effect, a certain (weak) sensitivity of 
L 2 on the sign of U, and especially the value of a, indi- 
cate that the approximations employed in Jl],|| were not 
accurate enough. We note that for the smallest disor- 
der values considered (W = 0.7,0.8) Li is of the order 
of the system size. Nevertheless a strong enhancement 
effect prevails. This interesting observation supports the 
claim [^) that persistent currents (in metallic systems) 
are enhanced by interaction effects. 

Bag model with medium range interaction. To further 
substantiate our results, we turn our attention to the 
"bag model" already mentioned in In the infinite 
square lattice in which each point (ni,^) corresponds 
to a wave function ip ni ,n 2 of t ne two-electron system, we 
consider now the transfer along the diagonal instead of 
the ni-direction. We associate with each center-of-mass 
coordinate n — n\ + rii a "tranverse" vector <p n con- 
taining Nb two-particle configurations with fixed n, but 
increasing distance \ni — n®]. In contrast to the previous 
study (where we have ignored the symmetry of the wave 
function) we assume VVii,n 2 — l 4 , n 2 .n 1 and consider only 
points in the upper triangle of the scheme. The coupling 
between the 4> n can be written as 

E(f> 2n = h 2n (j)2n + (02n+l + 4>1n-l) j 
E(j>2n+1 = /l2ri+102n+l + (02n+2 + 4>n) ■ (6) 



The form of the matrices hi and £( e >°) follows directly 
from (|l|) and we have to distinguish between fan and 
4>2n+i since only the former has an entry on the diagonal 
rii = n 2- The corresponding 2AT& x 2^ transfer matrices 
are given by 



M; 



(e,o) _ ( t^-^E-hi) 
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They have to be iterated on a set of Nb different start vec- 
tors — {4>u \ <t>d) w hich have to be orthonormalized 
from time to time as the iteration progresses Sum- 
ming up the logarithmic increases of the norms of these 
vectors one calculates the set of A^ Lyapunov exponents 
ji = of this transfer matrix. The start vectors 
tend to the stationary distribution of the eigenvectors of 
M T M, where M = JJl 



Mi e) Mi o) . 
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Unfortunately, the truncation of the transverse vectors 
(pi in the bag model introduces an artificial, long-range 
interaction between the two electrons: They feel a hard 
wall whenever their relative separation An reaches its 
maximum value (An) mnl = 2Nb- A detailed analysis of 
the bag model therefore meets with two problems: (i) to 
determine the localization length L\ for the noninteract- 
ing problem and (ii) to attribute any observed enhance- 
ment of L2 = £1 > Li either to the (unphysical) bag 
interaction or to the (physical) Hubbard interaction. 

We have performed numerical calculations for E = 
and N b = 100,200,300 in the disorder range 0.7 < 
W < 2.0. In our model the two electrons interact with 
[7 = 1 whenever their distance An is smaller than a fixed 
value (An) mi . This modification of the Hubbard inter- 
action is designed to enhance the physical effect as com- 
pared to the unphysical bag interaction. We have chosen 
(An) int = 20 (to be compared with (An) mnl = 2N b = 
200,400,600). In all cases the full Lyapunov spectrum 
and the corresponding modes (i.e. the average of the 
vectors <p\ over the full iteration process) have been cal- 
culated up to a statistical accuracy of at least 5 percent. 

Modes with dominant weight at very small An (max- 
imal An) owe their existence to the physical Hubbard 
interaction (the bag interaction). Therefore Fig. 3, where 
we show the modes corresponding to the two smallest 
Lyapunov exponents 71 and 72, demonstrates that the 
medium-range Hubbard-type interaction is slightly more 
effective in assisting coherent transport than the artifi- 
cial bag interaction. The dependence of L2 on the disor- 
der parameter W is determined, as for the previous case 
of finite systems, by plotting log £2 = l°g(l/7i) versus 
logVF. The corresponding curve for Nb = 300 is shown 
in Fig.2. We find L 2 cx W~ 2a with 2a = 3.42. This is in 
relatively good agreement with the previous model and 
confirms that L2 does not scale as L\ . 

Finally, we compare in the inset of Fig. 3 the small- 
est 36 Lyapunov exponents (for Nb = 100) with and 
without the Hubbard-type interaction. Two points de- 
serve to be mentioned: (i) The interaction selectively re- 
duces the smallest Lyapunov exponents, demonstrating 
the profound effect on the localization properties of the 
system, (ii) The smallest exponents which are no longer 
appreciably influenced by the interaction can serve (via 
L = 1 /j) as a reasonably good estimate for the localiza- 
tion length L± of the noninteracting system. The corre- 
sponding modes are the first ones to propagate without 
taking advantage of interaction effects. 

Q-statistics and scaling. To understand the discrep- 
ancy between our calculated and previously estimated 
exponents, we return to the overlap matrix element (|^). 
We have calculated the Q-distribution (see Fig. 4) for a Id 
system with W = 1.0 and N = 200 in the following way. 
Defining two adjacent boxes of size L\ ss 25 /W 2 = 25 
in the middle of the system, we selected all (noninteract- 
ing) eigenstates having their maximum in either of these 



blocks. Then the distribution was determined by con- 
sidering all Qij^im involving two wavefunctions from each 
block, respectively. We find deviations from a Gaussian 
behavior in two respects. First, while the off-diagonal 
terms Qijj m (i 7^ j and/or I 7^ m) are distributed sym- 
metrically around zero, the diagonal elements Qujj are 
always positive. These latter contributions describe si- 
multaneous transitions of two electrons (with opposite 
spins) between state i and state j. Therefore, they in- 
volve only two wavefunctions (and not four), are typi- 
cally larger than the off-diagonal terms and preserve the 
sign of the interaction U. Second, both the symmetric 
off-diagonal and the asymmetric diagonal parts of the 
Q-distribution are very sharply peaked at zero and have 
small but very long tails. Fig. 4 suggests that the par- 
ticular role of the diagonal elements and the strongly 
non-Gaussian character of the distributions might be re- 
sponsible for the modified exponent a. The form of p{Q) 
affects both the arguments given by Shepelyansky and 
by Imry. The band matrix model assumes that all pairs 
of two-particle states localized within a distance L\ are 
coupled with uncorrelated, uniformly distributed inter- 
action matrix elements. The scaling block picture esti- 
mates the interblock coupling using essentially the same 
assumption. In reality, sizable coupling matrix elements 
are sparsely distributed (and most likely quite correlated) 
within L\, so that the notion of a typical interblock cou- 
pling (or a well-defined bandwidth in the band matrix 
model) is not unproblematic. 

We nevertheless make an attempt to improve previous 
estimates. Following Imry @ we consider two blocks of 
size Li, each containing L\ two-particle states. We want 
to use the actual distribution p(Q) and therefore have to 
reduce a variable with a highly non-trivial distribution 
to a "typical" value needed for the Thouless scaling block 
picture. For illustration, we consider two extreme cases. 
First we calculate a restricted average (Q 2 ) r es of the over- 
lap matrix elements Q coupling a state in the first block 
and its nearest (in energy) neighbour in the second block. 
Such an average is unaffected by the particular distribu- 
tion of the Qujj, which couple levels with spacings of 
the order of 1/Li (instead of 1/L 2 ), and turns out to be 
consistent with L2 oc L 2 . Second, we calculate an un- 
restricted average over the full spectrum, including also 
large energy separations. We obtain (Q 2 ) oc LJ" 3 ' 30 which 
eventually gives L2 oc L} 70 in good agreement with our 
independent transfer matrix studies. We believe that the 
second averaging procedure is more appropriate. For a 
noninteracting system divided into blocks of size L\ the 
Thouless energy associated with a block is of the order of 
the one-particle level spacing <5i . The level uncertainty is 
therefore large compared to the two-particle spacing 82 
and should suffice to justify an unrestricted average B. 

In conclusion, we have found that the pair localiza- 
tion length L2 scales with an exponent a ~ 1.65, for 
a large range of system parameters and in contrast to 
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previous predictions. The deviation from the expected 
value could be traced back to the particular distribution 
of the overlap matrix elements Q, composed of a sym- 
metric off-diagonal and a positive diagonal part, both 
sharply peaked with long tails. By means of Imry's scal- 
ing picture the anomalous exponent, characterizing the 
interacting electrons, could be connected to properties of 
the noninteracting system. The relevance of the Qujj 
raises some doubts concerning the independence of the 
effect on the sign of the interaction. Studying the Q— 
distribution for d > 1 might be a useful tool to investi- 
gate and understand the effect discovered by Shepelyan- 
sky also in higher-dimensional systems. 
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FIG. 4. The distribution of overlap matrix elements Q for 
(i) a Id chain with N = 200, W = 1.0 (solid: diagonal terms, 
dashed: off-digonal terms) and (ii) a normalized Gaussian 
(dotted) with the same variance as (i). Inset: the same on a 
logarithmic scale. 
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FIG. 1. The five smallest and five largest inverse decay 
lengths ji(N) for a given sample with N = 100, W = 1, and 
E = 2 as a function of the Hubbard interaction U. 



FIG. 2. Log-log plot of the pair localization length L2 
as a function of the disorder W for (i) the finite chain with 
N = 100 (solid) and (ii) the bag model with N b = 300 (dot- 
ted). Both curves were calculated for (7 = 1 and E = 0. The 
exponent 2a — 3.3 and the behavior of L\ are also indicated. 



FIG. 3. The modes <j>ii corresponding to the two smallest 
Lyapunov exponents 71 (solid) and 72 (dashed) for the bag 
model with N b = 200, W = 1.5, U = 1, and E = 0. Inset: 
The smallest 36 Lyapunov exponents for the same model, but 
N b = 100 and (i) U = 1 (solid) and (ii) U = (dotted). 
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